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COMPUTER PROGRAM FOR THE CALCULATION OF MULTICOMPONENT 

CONVECTIVE DIFFUSION DEPOSITION RATES FROM CHEMICALLY 

FROZEN BOUNDARY LAYER THEORY* 

Suleyman A. Gokoglu* 

Analex Corporation 
Cleveland, Ohio 44135 

and 

Bor-Kuan Chen** and Daniel E. Rosner*** 

Yale University 

Chemical Engineering Department 
New Haven, Connecticut 06520 


SUMMARY 

The computer program developed based on multicomponent chemically frozen 
boundary layer (CFBL) theory for calculating vapor and/or small particle de- 
position rates Is documented. A specific application to perimeter-averaged 
Na 2 S 04 deposition rate calculations on a cylindrical collector Is demon- 
strated. The manual Includes a typical program Input and output for users. 


I. INTRODUCTION 

A wide variety of different engineering flelds-lncludlng chemical coating 
of metals, solid state electronics device fabrication, filtration, aircraft 
Icing, corrosion and fouling of gas turbine blades, heat exchangers, etc. -deal 
with vapor and particle capture, either to exploit or to suppress It, depending 
on the context. A comprehensive but tractable convective diffusion deposition 
rate theory has been developed at Yale University based on the assumption of a 
multicomponent "chemically frozen" boundary layer (CFBL) as outlined In Rosner 
et a_h [1]. The theory which makes full use of available transport coefficient 
and property Information Is not only applicable under conditions of multi- 
component vapor transport (e.g., CVD applications), but It Is also Intrinsi- 
cally capable of simultaneously dealing with particle transport provided the 
particles are small enough to be considered heavy molecules. The purpose of 
this manual Is to document a simple but useful computer code based on CFBL 
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theory for calculating alkali sulfate deposition rates. We hope that the 
manual will prevent further misapplications of CFBL theory as has been en- 
countered In the literature [ 2 ]. Detailed description of this program, with 
specific application to perimeter-averaged Na2S04 deposition rate calculations 
on a cylindrical collector Is given below. Typical program Input and output 
Is Included with explanations. The manual contains the program listing In 
Appendix C. 


SYMBOLS 

B7 thermophoretlc parameter, equation (35) 

C nozzle discharge coefficient 

C m h coefficient defining relative sensitivity of mass transfer and heat 

transfer to variable properties, equation (37) 

Cp heat capacity per unit mass at constant pressure 

c v heat capacity per unit mass at constant volume 

D Pick diffusion coefficient 

d collector (target) diameter 

dj nozzle jet diameter 

F function 

f fuel to air mass flow rate ratio 

I mainstream turbulence Intensity 

K coefficients of curve-fit In equation (11) 

k Boltzmann constant 

L mainstream turbulence length scale 

Le Lewis number (ratio of Pick to thermal dlffuslvlty) 

M molecular weight 

M average molecular weight of the mixture 

m" deposition mass flux 

Nu Nusselt number 

Nu Perimeter-averaged Nusselt number 

P pressure 

Pr Prandtl number (ratio of kinematic viscosity to thermal dlffuslvlty) 

R universal gas constant 

Re Reynolds number 

S spin parameter, equation ( 33 ) 

s defined by equation ( 31 ) 

Sc Schmidt number (ratio of kinematic viscosity to Flck dlffuslvlty) 

T absolute temperature 

U velocity 

v sonic velocity 

W a air mass flow rate 

X mole fraction 

Gre ek Symbols : 

<*7 thermal diffusion factor 

a T,oa thermal diffusion factor at Infinite temperature 

a 7 ,-l coefficient of T -1 In power series expansion of aj 

Lennard-Jones molecular Interaction energy (well-depth) parameter 
specific heat ratio, c p /c v 
thermal conductivity 
v viscosity 

v number of Na atoms In species 1 


2 


$ 

p 

a 


Wllke-Wasslljewa coefficients, equation (21) 
density 

Lennard- Jones molecular size parameter 
Gi angular speed of the collector (target) 

Q collision Integral 

Subscripts : 

av averaged quantity 

cl centerline 

cp constant properties 

D pertaining to FI ck diffusion coefficient 

h heat transfer 

1 species 1 

j jet exit plane, location "j" In figure 1 

k species k 

m mass transfer 

mix gas mixture 

o stagnation conditions 

t mainstream turbulence 

w target surface (wall) 

V pertaining to viscosity 

« mainstream, location "®" In figure 1 

Miscellaneous : 


CEC chemical equll Ibrlum code 

CFBL chemically frozen boundary layer 

CVD chemical vapor deposition 


II. DESCRIPTION OF PROGRAM WITH APPLICATION TO Na 2 S0 4 DEPOSITION 

The CFBL computer code was developed for and Is being used to Interpret 
Na 2 S0 4 deposition experiments at NASA Lewis Research Center, Cleveland, Ohio 
[3] because of Its Importance In hot corrosion. Therefore, here we specifi- 
cally concentrate on Na 2 S0 4 deposition rate prediction via the transport of 
many Na containing species to the collector, although the theory Is general In 
that other salts can easily be Incorporated. For a typical burner rig and 
collector configuration, and station nomenclature shown In figure 1, where any 
Na containing salt Is seeded Into the combustor, the deposition rate of Na 2 S0 4 
can be shown to be: 


m, 


Na 0 SO. 
2 4 


Na 2 S0 4 


F(turb) 


p ® 

d • 2^ 


l i..- N Vi * F i (Soret) 


(x xix b t.i / Le i.w \° 6 

i ,°° ' A 1,w> " A 1,w • F^Soret) * ^Le, J 



( 1 ) 


The multicomponent nature of the problem Is reflected by the appearance 
of summations which Include each Na carrier species 1. The present program 



considers four different possible Na carriers: Na, NaOH, Na 2 S0 4 and NaCl. 

Note the Implicit assumption ( checked below via the element Flux ratio con- 
tralnt) that once Na arrives at the surface of the collector by any of Na 
carriers, there Is always enough sulfur available (e.g., supplied to the system 
by fuel) to form Na 2 S 04 . 

The mole fraction of each Na carrier species In the free stream and at the 
collector surface (X< „ and X^ >w ) Is supplied to the program by the user. 

The theory leads to simple algebraic flux relations for each species since It 
assumes that there Is thermo-chemical equilibrium at stations and "w", 
although the boundary layer Itself Is chemically frozen. In the applications 
to Na 2 S 0 4 deposition experiments at NASA LeRC, the Chemical Equilibrium Code 
(CEC) developed at NASA LeRC [4] Is used to obtain the equilibrium compositions 
after combustion at specified temperatures and pressures. However, the CFBL 
code Is by no means restricted to coupling with CEC, and the user may provide 
Xi and X^ jW values by any other means available. 

Conditions prevailing at station "w" require the Imposition of one addi- 
tional constraint other than vapor/condensate equilibrium [1,5]. Because 
Na 2 S 04 Is being deposited, the molar fluxes of Na and S must stand In the 
ratio 2:1 at the prevailing vapor compositions. The element mass fractions of 
the Inorganic constituents characterizing the free stream are dictated by the 
experimental conditions; however, the situation at station "w" Is more compli- 
cated because transport across the boundary layer contributes to element 
"segregation." Therefore, the element mass fractions at the surface of the 
collector generally will not be the same as those In the free stream [ 5 ]. 

Hence, the CFBL program also calculates the total S molar flux to the surface 
with the same formulation given In equation (1). The dominant S carriers to 
the surface are chosen to be S0 2 , SO 3 , and H 2 S. Therefore, the user has to 
provide X ^ >00 and X^ w values for S0 2 , SO 3 , and H 2 S as Inputs to the pro- 
gram. For any particular set of X ^ >00 and X^ w values, obtained from equi- 
librium calculations and Input to the CFBL program, the output gives the cal- 
culated Na to S molar flux ratio. In order to satisfy both the thermodynamic 
(chemical equilibrium) and the transport (molar flux ratio) constraints, the 
user has to Iterate on X^ >w values obtained from equilibrium calculations 
with clever choices of element mass fractions. Finally, the properly chosen 
element mass fractions will give the proper equilibrium compositions of Na and 
S carriers at station "w" which In turn will satisfy the 2:1 Na to S molar 
flux ratio constraint. Experience with the CFBL program Indicates that In 
cases where the system contains excess S, the predicted Na 2 S0 4 deposition 
rates are practically Insenslstlve to the fulfillment of the transport re- 
quirement except close to the dew point temperature. + 

Modification of the CFBL program for application to the multicomponent 
chemical vapor deposition of condensates other than Na 2 S0 4 Is conceptually 
straightforward. It requires the proper choice of dominant carrier species and 
provision of relevant thermodynamic and transport parameters to the program. 
Generalized versions of the present CFBL code are being developed and applied 
to Na 2 S0 4 + K 2 S 0 4 and Na 2 S0 4 + K 2 S0 4 + L1 2 S0 4 solution condensate deposition 
situations used at Yale University [5] and NASA LeRC. 


tFor the discussion of the fact that "dew-point" Is not a pure thermodynamic 
quantity, the reader Is referred to references [ 5 ] and [ 6 ]. 
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III. EVALUATION OF THERMODYNAMIC AND TRANSPORT PROPERTIES 

Description of the property calculations given below assumes that the 
following Information Is supplied to the program by the user as Input for a 
typical configuration given In figure 1: 

dj jet exit diameter, cm 

d diameter of cylindrical collector (target), cm 

f fuel to air mass flow rate ratio 

Pj jet exit (ambient) pressure, atm 

P 0 stagnation pressure, atm 

T 0 stagnation temperature, K 

T w collector surface (wall) temperature, K 

W a air mass flow rate, g/sec 

Note that we are assuming a circular jet and a cylindrical collector. For 
other geometries proper characteristic dimensions must be supplied to calculate 
the jet exit area, the Reynolds number and mass transfer Nusselt number. 

The combustion Is assumed to be complete and dominated by the reaction: 

CH 2 + 3/2 0 2 -* C0 2 + H 2 0 (2) 

Hence, starting with air of Initial composition: 79.05 mole percent No and 

20.95 mole percent 0 2 , we estimate the composition of combustion product 
gases by: 


x _ 0.7905 

N 2 = 1 + 1 .0331 f 


(3) 


„ 0.2095 - 3.0993f 

0 2 = 1 + 1 .0331 f 


(4) 


x x 2.0662f 

C0 2 = a H 2 0 = 1 + 1 . 0331 f 


(5) 


All other species present In the combustion product gas stream are assumed to 
be of trace amounts and are neglected In the property calculations. 


A. Calculation of M, y, Tj and Uj 
Average molecular weight of the combustion gases Is obtained from: 


M = x n 2 m n 2 + x o 2 h o 2 + x co 2 m co 2 + x h 2 o h h 2 o 


( 6 ) 


where X Is the mole fraction of each Indicated species. 

The ratio of heat capacity at constant pressure to heat capacity at con- 
stant volume of the mixture Is calculated by: 
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(7) 


C p.m1x _ c p.m1x 

Cv ' m1x c p,mix - < R/ "> 

where R Is the universal gas constant (1.9872 cal/mole - K), and the calcula- 
tion of Cp >rn ^ x is described below. Assuming an isentropic jet, the tempera- 
ture and the velocity at the jet exit plane is given by: 


T 


j 



GS 


(Y-1)/y 


( 8 ) 


and 




W a (l + f) 

(■?)•■ 


(9) 


where R is the universal gas constant (82.057 atm-cm 3 /mole-K) , and C is the 
discharge coefficient of the jet nozzle to be supplied by the user. The dis- 
charge coefficient depends on the geometric shape of the nozzle and the 
Reynolds number based on the nozzle diameter. If no information is available 
for the nozzle discharge coefficient, then the program automatically assigns 
unity as the value of C. Note that the velocity calculated from equation (9) 
based on the total mass flow rate and the effective area of the nozzle exit 
should agree with the velocity which can be calculated also from the prevail- 
ing Mach number based on the isentropic jet assumption as: 



(9a) 


where v is the sonic speed in the gas mixture at T 0 . In equation (9) it is 
assumed that the mixture obeys the ideal gas law, so that the density of the 
mixture is simply: 


p 


rt j 


( 10 ) 


B. Calculation of c p , x, y, and 0^ 

The dimensionless heat capacity of each dominant species is computed from 
a curve fit: 



K 3.1 tZ * K 2., T * K 1.1 


(11) 


where R is the universal gas constant (1.9872 cal/mole - K). Two different 
sets of constants, K, are used for T > 1000 K and T < 1000 K for a better fit. 
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Viscosity, v], and thermal conductivity, of dominant species are cal- 
culated using Chapman-Enskog theory: 


^ = 2.6693x10 



x - — 
*1 " H 1 


i5 >K32 (y!i.|) 


(g/cm-sec) 


p^ (cal/cm-sec-K) 


( 12 ) 


(13) 


where Is the Lennard-Jones molecular size parameter of the species and 
fl»,1 Is the colllson Integral for viscosity which Is a function of the 
dimensionless temperature kT/e^ where k Is the Boltzmann constant and 
the Lennard-Jones molecular Interaction energy (well-depth) parameter for 
species 1. 

lhe binary FI ck diffusion coefficient of trace species k In each of 
dominant species 1 Is also calculated from Chapman-Enskog theory: 


D kJ = 0.0018583 


/Wk) 


Po k1 fl D,k1 


(cm /sec) 


(14) 


where 


°ki = 2 * (o k + °1 } (A) (15) 

and Is the collision Integral for diffusion coefficient which (just 

like «’j) Is also a function of kT/c^^ where: 


fki tSh 

k = k 


(16) 


The dimensionless thermal diffusion factor for trace species 1 Is convenient- 
ly expressed below after a curve-fit to calculations based on Chapman-Enskog 
theory. 


T , 1 = a T,«,1 


T.-1.1 


(17) 


The resulting least-squares constants, aj „ ^ and ay _ 1 * for a number of 
species are given In reference [7]. ’ ’ ’ ’ 


It should be noted that unless experimentally determined, the uncertain- 
ties Involved In the estimation of a and e/k are directly reflected In p, 
X, D, and aj values [8]. Therefore the accuracy of prediction of the CFBL 
program Is limited by the accuracy of such transport parameters. 
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C. Calculation of Mixture Properties 
Specific heat of the gas mixture Is simply given by: 

c p,m1x = ^ IV. * * * * * X 1 C p,1 (18) 

where Is the mole fraction of species 1. Viscosity and thermal conduc- 
tivity of the mixture are calculated from similar formulae 


v m1x 


Nnlx 



1 



Vi 

? Xk$k1 


Vi 


? X k*k1 


(19) 


( 20 ) 


where §|^ are the Wllke-Wasslljewa coefficients given by: 


§ 


kl 




( 21 ) 


The diffusion coefficient of species 1 In the mixture Is calculated from: 


°1 .mix 


The reader Is referred to reference 
transport property calculations. 



( 22 ) 


[9] for the details of thermodynamic and 


IV. DIMENSIONLESS NUMBERS 

The dimensionless numbers given below are evaluated either at station "®" 

or "w" depending on relevance. The static temperature at station Tc, Is 

assumed to be equal to Tj, calculated In Section III. A. The collector surface 

temperature (or a representative average temperature If the collector Is not 

Isothermal) and the ambient pressure, Pa,(=Pj), are supplied as Inputs by the 

user. 
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(23) 


Pr 


c p.m1x y m1x 

Nnlx 


- y m1x 
Sc 1 " ~ 


1 .mix 


(24) 


Le, = 


Pr 

Sc 


1 .mix 


i 


[x m1x /(pC p,m1x )] 


(25) 


pU d 

Re = — — (26) 

v m1x 

Calculation of the Reynolds number requires knowledge of U^. Although 
Uj Is obtained In Section III. A, U^ Is not necessarily the same as Uj de- 
pending on the jet velocity profile and the jet divergence angle. For jet 
diameters larger than the collector diameter (usually) the case), one needs the 
centerline velocity at station "®" for the proper calculation of Re. However, 
Uj calculated In Section III. A. Is the average velocity at station "J". 
Therefore, two corrections are needed; first, to correct average Uj to center- 
line Uj, F (shape), and second, to correct centerline Uj to centerline U,,,, 

F ( d 1 v ) . These corrections take Into account the boundary layer effects and 
the flow area variation from station "j" to station depending on jet diver- 
gence angle. Hence: 


U^ = Uj • F (shape) • F (dlv) 


(27) 


where F (shape) and F (dlv) are defined as: 


F (shape) = 



(28) 


F (dlv) = li^ (29) 

cL,j 

F (shape) and F (dlv) are optional Inputs to the program. If no Information 
Is available, the program automatically assigns unity for both. Based on 
available data for fully developed turbulent flow In smooth circular ducts, 
a suggested estimate for F (shape) Is: 

F (shape) s ( s + U (2s + 1) ( 30 ) 

2s* 


where 
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(31) 


s 



The mass transfer Nusselt number correlations depend, among other things, on 
the geometry of the collector used. In the NASA LeRC experiment [3], a cylin- 
drical collector was used. The example program listing given at the end of 
this manual also assumes a cylindrical collector. By making the replacement 
Pr -» Sc^, the perimeter averaged mass transfer coefficient, Nu m *, becomes 
based on the heat transfer correlations recommended for stationary Inflnltly 
long cylinders [10]. The additional effects of variable properties and Mach 
number are approximated by Incorporating a post-multiplier term based on the 
recommendations of references [11] and [12]. 



(0.40 Re 1/2 + 0.06 Re 273 ) 



(32) 


It should be noted that If a different collector geometry Is used and/or the 
surface roughness of the collector can have an appreciable effect on mass 
transfer, the equation (32) may not be applicable, and the user has to modify 
the program to supply the best available corresponding mass transfer Nusselt 
number correlation from the literature. 


In many experiments, uniformity of collector surface temperature Is ob- 
tained by collector rotation and this requires additional modifications of the 
formulation given by equation (1). In the NASA LeRC deposition experiments the 
rotational speeds of the smooth collectors used have been neglected based on 
reference [13], which suggests that If one defines a "spin" parameter 


_ ud 
“ 211 


(33) 


where o Is the angular speed of the target, then equation (1) (or eq. (32)) 
should be corrected by the multiplier 


„ 0.315 

F (spin) s (1 + 2S ; 


(34) 


V. THERMAL (SORET) DIFFUSION EFFECT, F (SORET) 

Inclusion of thermal (Soret) diffusion In the mass transfer calculations 
Introduces a nondlmenslonal thermophoretlc parameter, [14], which can be ex- 
pressed as: 



-a 


T, 1 ,w 



(35) 


The effects of thermal (Soret) diffusion ("thermophoresis" for small particles) 
Is thoroughly discussed and correlated In references [7, 14-17] for both vapor 
and small particle mass transfer. For the current applications of the CFBL 
program to NASA LeRC deposition experiments, the emphasis has been on vapor 
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deposition. Therefore, the "sink" effect of thermal diffusion Is justifiably 
neglected, and only the "suction" effect Is considered, l.e.: 


F 1 (Soret) £ F 1 (suction) £ 1 _ eX p*j B j (36) 

T* 1 

For small particle mass transfer applications where Inertial effects are not 
Important, the user should modify the program to Include the "thermophoretlc 
sink" effect as described In references [14] and [16]. For such applications 
a method for estimating small particle transport properties (D^, aj etc.) 

Is given In reference [17]. 


In the program the user also has the option of completely "turning off," 
If so desired, the thermal diffusion effect from Inclusion In the calculations 
by specifying the Input variable "FSORET" as "false". In this simple way we 
can demonstrate the relative Importance of thermal diffusion mass transport In 
any multicomponent CVO situation. 


It Is known that variable properties across the boundary layer affect 
mass transfer rates In a somewhat different manner than they do the better 
known heat transfer rates. In case of thermal diffusion, this effect Is in 
addition to the one already Incorporated In the perimeter-averaged Nusselt 
number (eq. (32)). Certainly, the effect will be different depending on the 
temperature dependence of the transport properties of the species In question 
and the characteristics of the mass transfer boundary layer. Therefore, It Is 
convenient to Introduce a new dimensionless parameter, C m ^ j, defined as: 


(Nu /Nu ) 
' m m,cp' 


1 


'mh,1 ~ (Nu./Nu. ) 
h h,cp 


(37) 


Simple engineering correlations to predict the effect of variable properties on 
mass transfer (the numerator In eq. (37)) and on heat transfer (the denomenator 
In eq. (37)) are given In reference [12]. However, If the space reserved for 
c mh,1 the Input data Is left blank, the program automatically assigns a 
value of unity for each C m h 


VI. MAINSTREAM TURBULENCE EFFECT, F (turb) 

It Is known that mainstream turbulence, apart from turbulence within the 
boundary layer, Influences the transport rates to surfaces [18-19]. Although 
there Is considerable quantitative disagreement among various researchers, the 
CFBL program uses the correlation given In reference [20] based on perimeter- 
averaged heat transfer measurements for cylinders In cross flow. If one has 
knowledge of free stream turbulence Intensity, 1^, and axial macroscale of 
turbulence, L^, then the correlation at a particular Re for a cylinder of 
diameter d Is given In the separable form: 

F (turb) £ 1 + fct^^Re) • fct 2 (L t /d) (38) 

where f c t-j and fct 2 are presented In graphical form In reference [20]. The 
CFBL program uses the following curve fits for f ct-j and fct 2 * 
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fct-| ( 1 1 Re ) 


fct 2 ( L t /d) 


r 

I 

V e \ 3/2 1 

12.375 

1 M 

- -V) 


\ 

io 4 / 

{ 




If I t Re < 10 4 


9.0 + 3.375X10" 4 I t Re If I t Re > 10 4 

L t /d > 2 
L t /d < 2 


0.124xl0" 3 [(L t /d) - 11.0] 2 + 2.0X10" 3 If 
-4.0xl0" 3 [(L t /d) - 1.75] 2 + 12.25X10 -3 If 


(39) 


(40) 


F (turb) probably also depends on Sc 1t but the Information Is not yet suffi- 
cient In the literature to Incorporate It Into the present formalism. The user 
has three options In using the F (turb) factor: 


1. If one chooses not to use the F (turb) factor, then leave the variable 

"TURB" blank In Input data. In which case the program automatically 
assigns the value of unity for F (turb). 

2. If the particular F (turb) to be used by the program Is known, then 

that value should be read In as Input for the variable "TURB", but 
leave the variables "TURIN" and "TURL" blank. 

3. If one has knowledge of Ij- and and chooses to use the correla- 
tion supplied by the program based on equations (38) to (40), then 

Input data for "TURB" Is dummy, but "TURIN" (percent 1+ = I+/100) and 
"TURL" Input values must be read In. 


below 

MAIN 

PARAH 

CALESL 

CALMIX 

COLINT 


MIXPRO 

C0E0IF 

DIF 

TMDF 

NUM 

TURBL 


VII. ROUTINES 

The names and functions of the subroutines used In the program are given 


calculates the deposition rate based on equation (1) 
provides e/k, a, and molecular weight of dominant species In the 
gas mixture 

provides the dimensionless heat capacity of dominant species In the 
gas mixture from a curve fit 
calculates c p>m1x In cal/g-K 

evaluates ’and Qq from a curve fit to the already supplied 
kT/c versus and fl D data [8]. Also calculates and x< 
for dominant species 
calculates ii m ^ x and x,^ 

0(^ of trace species k In 1 

D 1 .mix 
F (Soret) 

Nu m,1 

F ( turb) 


calculates 

calculates 

calculates 

calculates 

calculates 
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VIII. DESCRIPTION OF PROGRAM INPUT 


The Input data to the CFBL program Is read In three different stages. 

The first stage reads In one card which contains the following Information 
according to the format specified. 

READ( 5, 1 ) IRUN, ITYPE, F, WA, TW, TO, PO, PJ, DC, WOBS, FSORET 
1 FORMAT ( 212 , 8F9.0, 3X, LI) 

IRUN The user can use this variable for his own purposes (e.g., to keep 
track of which experiment, which run, etc.). 

ITYPE lhe user can use this variable for his own purposes (e.g., to keep 
track of the type of salt fed to the combustor, etc.). 

F Fuel to air mass flow rate ratio 

WA Air flow rate, g/sec 

TW Collector (target) surface temperature, K. Give an average tempera- 

ture If the collector surface Is not Isothermal. 

TO Stagnation temperature, K 

PO Stagnation pressure, atm 

PJ Jet exit (ambient) pressure, atm 

DC Discharge coefficient of jet nozzle. Optional Input. Leave blank 

If no Information Is available. See Section III. A. 

WOBS Experimentally observed deposition rate. Used In calculating the 
percent error In the prediction, mg/hr. 

FS0RE1 Logical variable. Specify as 1 T 1 If thermal diffusion Is desired to 
be Included In the calculations, or as ' F 1 If otherwise. 

The second stage reads In one card which contains the following Information 
according to the format specified: 

READ( 5,998) DIAW, LW, DJ, TURB, TURIN, TURL, SHAPE, DIV 
998 FORMAT(8F10.0) 

DIAW Diameter of the cylindrical target, cm 
LW Length (height) of the cylindrical target, cm 

DJ Diameter of circular jet nozzle exit plane, cm 

TURB Turbulence enhancement factor, F (turb). Optional Input. 

(a) Leave blank If no consideration Is required. 

(b) Give your own estimate If Information Is available, 
but then leave "TURIN" and "TURL" blank, 

TURIN Percent mainstream turbulence Intensity (1-fVlOO). Optional Input. 

(a) Leave blank If no consideration Is required. 

(b) Leave blank If you supply your own "TURB". 

(c) Give your estimate If program consideration Is desired, 
but then specify "TURL", too. 

TURL Axial macroscale of mainstream turbulence, cm. Optional Input. 

(a) Leave blank If no consideration Is required. 

(b) Leave blank If you supply your own "TURB". 

(c) Give your own estimate If program consideration Is desired, 
but then specify "TURIN", too. 

SHAPE F (shape). Optional Input. Leave blank If no Information Is 
available. See Section IV. 

DIV F ( d 1 v) . Optional Input. Leave blank If no Information Is 

available. See Section IV. 
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The third stage reads in seven cards, one for each species. The order Is 
given by: 

1 = NaOH 

2 = Na 

3 = Na 2 S0 4 

4 = S0 2 

5 = S0 3 

6 = NaCI 

7 = H 2 S 

Note that this Input Is specific to Na 2 S0 4 deposition only, and must be modi- 
fied for other salts and/or salt solutions accordingly. Each card contains the 
following Information according to the format specified. 

READ( 5,2) XW(I), XJ(I), CMH(I) 

2 F0RMAT(2E15.4, F10.0) 

XW(I) Equilibrium mole fraction of species 1 at station "w" 

XJ(I) Equilibrium mole fraction of species 1 at station "j" or equiva- 

lently at station "®" 

CMH(I) Ratio of effect of variable properties on mass transfer to heat 
transfer. Optional Input. See Section V. 

A typical Input data set Is given In Appendix A. 


IX. DESCRIPTION OF PROGRAM OUTPUT 


The program output contains all the Information read Into the program as 
Input. Description of the additional variables printed out Is given below: 


D(I) 

NU( I ) 

SC(I) 

TAU(I) 
F(SORET) (I) 
M(I) 

YNAJ 

YNAW 

YSJ 

YSW 

AM 

TURB 

RE 

SUM 

SMS 

RATIO 

ERROR 

GM 

PR 


Diffusion coefficient of species 1 In gas mixture, D^ m ^ x . 

In cm^/sec * 

Mass transfer Nusselt number of species 1, Nu m ^ 

Schmidt number of species 1, Sc 1t at station "i" 

Thermophoretlc parameter of species 1, Bt * 

F^(Soret) of species 1 ’ 

Equivalent Na 2 S0 4 deposition rate of species 1 mass transfer 
rate In mg/hr 

Total Na mole fraction In flame 
Total Na mole fraction at station "w" 

Total S mole fraction In flame 
Total S mole fraction at station "w" 

Average molecular weight of gas mixture (neglects trace species) 
F(turb) 

Reynolds number at station "•»“ 

Total Na mass flux to the surface In g/cm^-sec 
Total S mass flux to the surface In g/cm^-sec 
Na to S molar flux ratio at the surface 
Percent error between experimentally observed and predicted 
Na 2 S0 4 deposition rates 

Y = c p,m1x/ c v 

Prandtl number at station "®" 
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X(I) 

TJ 

RHOJ 

UJ 

E1AMIX 
LAM IX 
CMIX 


Mole fraction of dominant species In flame 
Jet exit static temperature, Tj = T*, K 
Density of gas mixture at station "j" or equivalently "« M In 
g/cnr 

Jet exit velocity, Uj, In cm/sec 
p m ^ x at station "®" In g/cm-sec 

at station "®" In cal/cm-sec-K 
c p,m1x at station "®" In cal/g-sec 


A typical 


output data set Is given In Appendix B. 
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123 4 56789012345678901234567890123456789012345678901234567890123 4 5678901234567890X XXX REFERENCE LINE 


. F WA 

19 2 0.0481 20.0214 

DIAW LW 

1.905 1.2649 

XWCI) 

6 . 981 0 E-14 
3.1280 E-20 
3 .6060E-12 
8.0742E-06 
1 .3345E-05 
0.0 

5.9210 E-31 


TW 

900 . 0 
DJ 

2.54 

XJ(I) 

6 . 38 1 0 E- 0 6 
8 .2740E-07 
2.6170E-13 
2.4994E-05 
4 .8236E-08 
0.0 

1 . 1600E-14 


TO 

1885.5 

TURB 

0.0 

CMH(I) 
0.0 
. 0 
. 0 
.0 
. 0 
0 
0 


PO 

1.068 

TURIN 

0.0 


PJ 

1.000 

TURL 

0.0 


DC W0B5 

1.000 27.00 
SHAPE 
0.0 


FSORET 

T 

DIV 

0.0 


123456 78 901234567890123456789012345678901234567890123456789012345678901234567890XXXX REFERENCE LINE 


APPENDIX A 


CFBL THEORY FOR SODIUM SULFATE DEPOSITION RATE 

ALL GAS PROPERTIES PERTAIN TO TJ,PJ 

IRUN =19 

ITYPE = 2 

F = 0.048100 

TU (K) = 900.000 

TO (K) = 1885.500 

PO (ATM) = 1.068000 

PJ (ATM) = 1.000000 

DC = 1.000000 

WOBS (MG/HR) = 27.0000 

FSORET =T 

WA (G/SEC) = 20.02139 

DIAW (CM) = 1.90500 

LW (CM) = 1.26490 

DJ (CM) = 2.54000 

TURIN = 0.0000 PERCENT 

TURL (CM) = 0.00000 

SHAPE = 0.0000 
DIV = 0.0000 


I=l=NAOH , I=2=NA , I=3=NA2$04 , I=6=NACL 


I D( I ) 

1 0.29905E 01 

2 0.37881E 01 

3 0 . 18316 E 01 

6 0.25178E 01 


NU ( I ) 

0.78633E 02 
0.71538E 02 
0 . 9567 0E 02 
0.84235E 02 


SCCI) 

0.10820E 01 
0.85418E 00 
0 . 17666 E 01 
0 . 12851 E 01 


X(I),J 
0. 63810 E- 05 
0.82740E-06 
0.26170 E-12 
0.00000 


XCI),W TAUCI) 

0.69810E-13 0. 73861 E- 01 

0.31280E-19 -0.23133E-01 

0.36060E-11 0.26800E 00 

0.00000 0.15291E 00 


FCSORET) , I 
1.03738 
0.98848 
1.13998 
1.07840 


YNAJ 

0.720840E-05 


YNAW 

0.728181E-11 


YSJ 

0 . 25 0422E-04 


YSW 

0.214192E-04 


AM = 28.81070 

TURB = 1.00000 

RE = 12913.625 

SUM (G/CM*X2/SEC) = 0.14063E-06 

SMS (G/CMXX2/SEC) = 0.26371E-06 

RATIO = 0.7437 

PREDICTED NA2S04 DEPOSITION RATE (MG/HR) = 11.83948 

ERROR (%) = -56.1501 

GM = 1.26601 

PR = 0.689647 

X(N2) = 0.753079 

X(02) = 0.057563 

X(H20) = 0.094679 

X(C02) = 0.094679 

TJ (K) = 1859.6157 


RH0J(G/CM**3) UJ(CM/S) ETAMIX( POISE) LAMIXCCAL/CM/K/S) CMIX(CAL/G/S ) 

0.188806E-03 0.219344E 05 0.610923E-03 0.290277E-03 0.327682E 00 


MCI) 

0.10364E 02 
0.14756E 01 
-0.64711E-05 
0.00000 


APPENDIX B 



APPENDIX C 


0000100 c 


0000200 C*X#X*XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX. 
0000300 C 
0000400 c 

0000500 C VARIABLE NAME LIST 

0000600 C 
0000700 C 

oooosoo c 

0000900 C VARIABLE DESCRIPTION 

0001000 c 

0001100 C ADT = TEMPORARY STORAGE VALUE FOR KT/E 

0001200 C A 1 1 , A J J = TEMPORARY STORAGE VALUE FOR DAT0S(3,I) 

0001300 C AJ = AREA OF JET NOZZLE (CM**2) 

0001400 C ALPAC79) = KT/E AS INPUT PERMANENT DATA 

0001500 C ALPHA ( I ) = THERMAL DIFFUSION COEFFICIENT FOR TRACE SPECIES I 

0001600 C AM = MEAN MOLEC. WEIGHT OF GAS MIXTURE! NEGLECTS TRACE SPECIES) 

0001700 C APHMIX = THERMAL DIFFUSIVITY OF MIXTURE AT STATION J 

0001800 C APHW = THERMAL DIFFUSIVITY OF MIXTURE AT STATION W 

0001900 C AIJ = AREA OF THE COLLECTOR ( TARGET ) , (CMXX2) 

0002000 C 

0002100 C C = CMIX INTERMEDIATE 

0002200 C CJ = DIMENSIONLESS HEAT CAPACITY 

0002300 C CMIX = MIXTURE HEAT CAPACITY AT STATION J (CAL/(G*K)) 

0002400 C CMHCI) = RATIO OF EFFECT OF VARIABLE PROP. ON MASS TO HEAT TRANS. 

0002500 C OPTIONAL INPUT. LEAVE BLANK IF NO INFORMATION AVAILABLE. 

0002600 C COEF1 = COEFFICIENTS FOR COMPUTING HEAT CAPACITY ( 1000K<T<5000K) 

0002700 C C0EF2 = COEFFICIENTS FOR COMPUTING HEAT CAPACITY ( 300K<T<1000K) 

0002800 C CONDT = THERMAL CONDUCTIVITY (TEMPORARY) , ( CAL/ ( CMXSECXK) ) 

0002900 C CP = MIXTURE HEAT CAPACITY AT STATION U (CAL/CGXK)) 

0003000 C 

0003100 CD = TEMPORARY STORAGE FOR BINARY DIFF. COEF. (CM**2/SEC) 

0003200 C DC = DISCHARGE COEFFICIENT OF JET NOZZLE. OPTIONAL INPUT. 

0003300 C LEAVE BLANK IF NO INFORMATION IS AVAILABLE, IN WHICH 

0003400 C CASE THE PROGRAM AUTOMATICALLY ASSIGNS UNITY AS ITS VALUE 

0003500 C DA = DIFF. COEFF. OF TRACE SPECIES AT STATION J ( CMXX2/SEC) 

0003600 C DATOS ( 1 , 1 ) = E/K FOR SPECIES I 

0003700 C DATOS ( 2 , I ) = SIGMA FOR SPECIES I (COLLISION DIAMETER PARAMETER) 

0003800 C DATOS (3* I ) = MOLECULAR WEIGHT OF SPECIES I 
0003900 C DEN = DENOMINATOR 

0004000 C DIAW = DIAMETER OF THE CYLINDERICAL COLLECTOR (CM) 

0004100 C DIFF = TEMPORARY STORAGE FOR BINARY DIFF. COEF. (CMXXZ/SEC) 

C004200 C DIV = DIVERGENCE FACTOR FOR NOZZLE JET DEFINED AS THE RATIO 

0004300 C OF CENTERLINE VELOCITY AT COLLECTOR LOCATION TO 

0004400 C CENTERLINE VELOCITY AT JET EXIT. DEPENDS ON THE JET 

0004500 C DIVERGENCE ANGLE. OPTIONAL INPUT, LEAVE BLANK IF NO 

0004600 C INFORMATION AVAILABLE. 

0004700 C DJ = DIAMETER OF JET NOZZLE EXIT PLANE (CM) 

0004800 C DW = DIFF. COEF. OF TRACE SPECIES AT STATION W (CMXX2/SEC) 

0004900 C 

0005000 C E/K = RATIO OF MOLECULAR INTERACTION ENERGY PARAMETER 

0005100 C (LENNARD-JONES) TO BOLTZMANN CONSTANT (K) 

0005200 C EOKI.EOKJ = TEMPORARY STORAGE FOR DATOS ( 1 , I ) =E/K 

0005300 C EOVKIJ = E/K, SQUARE ROOT OF EOKIXEOKJ 

0005400 C EPSOVK = DATOS ( 1 , I ) =E/K 
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0005500 C 
C 0 0 5 6 0 0 C 
0005700 C 
0005800 C 
0005900 C 
0006000 C 
0006100 C 
0006200 C 
0006300 C 
0Q06400 C 
0006500 C 
0006600 C 
0006700 C 
0006S00 C 
0006900 C 
0007000 C 
0007100 C 
0007200 C 
0007300 C 
0007400 C 
0007500 C 
0007600 C 
0007700 C 
0007800 C 
0007900 C 
OOOSOOO C 
0008100 C 
0008200 C 
0008300 C 
0003400 C 
0008500 C 
0003600 C 
0008700 C 
0003300 C 
0003900 C 
0009000 C 
0009100 C 
0009200 C 
0009300 C 
0009400 C 
0009500 C 
0009600 C 
0009700 C 
0009800 C 
0009900 C 
0010000 C 
0010100 C 
0010200 C 
0010300 C 
0010400 C 
0010500 C 
0010600 C 
0010700 C 
0010800 C 


ERN 

ERROR 

ET 

ETA 

ETAMJ 

ETAMIX 

F 

FSORET 


GM 

GMR 

G 

ITYPE 
I RUM 

l AMD A 

LAMP 

LAMIX 

LEJ 

LOD 

LW 

LEW 

M 

MCI) 


MU = 

OMEDT , OVDT = 
OMEDC79) = 

OME (79) = 

PHI = 

PJ = 

PO = 

PR = 

PROV = 

PSI = 

R = 

RATIO = 

RE = 

REJ = 

Rim j = 

R! !OW = 

SC = 

SHAPE = 


DIFFERENCE BETWEEN PREDICTED AND OBSERVED DEPO. RATES 
PERCENT ERROR BETWEEN PREDICTED AND OBSERVED DEP. RATE 
EXPONENTIAL (-TAU ) , SORET INTERMEDIATE 
MIXTURE VISCOSITY AT STATION W (G/(CM*SEC)) 

TEMPORARY VALUE FOR MIXTURE VISCOSITY (G/(CM*SEC)) 

MIXTURE VISCOSITY AT STATION J (G/(CM*SEO) 

FUEL TO AIR RATIO (BY MASS) 

LOGICAL VALUE PLACED INTO THE INPUT DATA FILE* IF IT IS 
SET TO *T f THEN THE SORET EFFECT WILL BE INCLUDED. IF IT 
IS SET TO *F* THEN THE SORET EFFECT WILL BE EXCLUDED. 

GAMMA 

(GAMMA-1 . )/GAMMA 
WASSIL J EWA -WILKE FACTOR 

NOT USED BY PROGRAM. OPTIONAL FOR USER PURPOSES. 

NOT USED BY PROGRAM. OPTIONAL FOR USER PURPOSES. 

MIXTURE THERMAL CONDUCTIVITY AT STATION W (CAL/C CM*SEC*K) > 
TEMPORARY THERMAL CONDUCTIVITY (CAL/(CM*SEC*K) ) 

MIXTURE THERMAL CONDUCTIVITY AT STATION J (CAL/(CM*SEC*K) ) 
LEWIS NUMBER FOR TRACE SPECIES AT STATION J 

RATIO OF MAINSTREAM TURBULENCE LENGTH SCALE TO COLLECTOR DIAM. 
LENGTH OF THE CYLINDERICAL COLLECTOR (CM) 

LEWIS NUMBER FOR TRACE SPECIES AT STATION W 

DAT0S(3,1)=M0LECULAR WEIGHT OF SPECIES I 

EQUIVALENT NA2S04 DEPOSITION RATE OF TRACE SPECIES MASS 

TRANSFER RATE (MG/HR) 

NUSSELT NUMBER FOR MASS TRANSFER 

COLLISION INTEGRALS 

COLLISION INTEGRAL FOR MASS DIFFUSIVITY 

COLLISION INTEGRAL FOR VISCOSITY OR THERMAL CONDUCTIVITY 

F(TURB) INTERMEDIATE RELATED TO TURBULENCE INTENSITY 
JET EXIT PLANE PRESSURE (ATM) 

STAGNATION PRESSURE (ATM) 

PRANDTL NUMBER AT STATION J 

INTERMEDIATE IN CALCULATION OF ETAMIX 

F(TURB) INTERMEDIATE RELATED TO TURBULENCE LENGTH SCALE 

GAS CONSTANTS. 9872 CAL/( GMOLE*K) 

SODIUM TO SULFUR MOLAR FLUX RATIO 
REYNOLDS NUMBER AT COLLECTOR LOCATION 
REYNOLDS NUMBER AT OUTLET OF JET 
DENSITY OF MIXTURE AT JET OUTLET 
DENSITY OF MIXTURE AT WALL 

SCHMIDT NUMBER FOR TRACE SPECIES AT STATION J 
SHAPE FACTOR (CENTERLINE VELOCITY)/( AVERAGE VELOCITY) 

OPTIONAL INPUT. LEAVE BLANK IF NO INFORMATION AVAILABLE. 
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0010900 

0011000 

0011100 

0011200 

0011300 

0011400 

0011500 

0011600 

0011700 

0011800 

0011900 

0012000 

0012100 

0012200 

0012300 

0012400 

0012500 

0012600 

0012700 

0012800 

0012900 

0013000 

0013100 

0013200 

0013300 

0013400 

0013500 

0013600 

0013700 

0013800 

0013900 

0014000 

0014100 

0014200 

0014300 

0014400 

0014500 

0014600 

0014700 

0014800 

0014900 

0015000 

0015100 

0015200 

0015300 

0015400 

0015500 

0015600 

0015700 

0015800 

0015900 

0016000 

0016100 

0016200 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SIGMA = MOLECULAR SIZE PARAMETER ( LENNARD-JONES) 

SIGI,SIGJ = TEMPORARY STORAGE FOR DATOSC 2 , I ) =SIGMA 
SIGK,SIGL = TEMPORARY STORAGE FOR DATOSC 2, I ) =SIGMA 
SIG2 = DAT0S(2,I)=SIGMA 

SIG = SIGMA VALUE FOR SPECIES I 

SMC I) = MASS FLUX OF TRACE SPECIES I (G/CMXX2/SEC) 

SMS = TOTAL SULFUR MASS FLUX TO THE SURFACE CG/CMXX2/SEC) 

SUM = TOTAL SODIUM MASS FLUX TO THE SURFACE (G/CMXX2/SEC) 

T = TEMPERATURE (K) 

TAU = SORET SUCTION PARAMETER FOR SPECIES I 

TO = STAGNATION TEMPERATURE (K) 

THERM = THERMAL DIFFUSION FACTOR FOR SPECIES I 
TJ = JET (STATIC) TEMPERATURE (K) 

TURB = TURBULENCE ENHANCEMENT FACTOR . OPTIONAL INPUT. 

A) LEAVE BLANK IF NO CONSIDERATION REQUIRED. 

B) GIVE YOUR OWN ESTIMATE IF INFORMATION AVAILABLE, BUT 
LEAVE ? TURIN 1 AND ' TURL ' BLANK IF YOU DO SO. 

TURIN = PERCENT MAINSTREAM INTENSITY. OPTIONAL INPUT. 

A) LEAVE BLANK IF NO CONSIDERATION REQUIRED. 

B) LEAVE BLANK IF YOU SUPPLY YOUR OWN f TURB ' . 

C) GIVE YOUR ESTIMATE IF PROGRAM CONSIDERATION DESIRED, 
BUT THEN GIVE AN ESTIMATE FOR 'TURL 1 , TOO. 

TURL = MAINSTREAM LENGTH SCALE (MACRO) , CM. OPTIONAL INPUT. 

A) LEAVE BLANK IF NO CONSIDERATION REQUIRED. 

B) LEAVE BLANK IF YOU SUPPLY YOUR OWN 'TURB'. 

C) GIVE YOUR ESTIMATE IF PROGRAM CONSIDERATION DESIRED, 
BUT THEN GIVE AN ESTIMATE FOR 'TURIN', TOO. 

TW = WALL TEMPERATURE (K) 

UJ = JET EXIT VELOCITY 

VISC(I) = VISCOSITY OF SPECIES (G/CCMXSEC)) 

V(I) = NUMBER OF SODIUM ATOMS IN TRACE SPECIES I 

WA = AIR FLOW RATE (G/SEC) 

WI,WJ = TEMPORARY STORAGE VALUES FOR DAT0S(3,I) 

WPRED = PREDICTED NA2S04 DEPOSITION RATE (MG/HR) 

WOBS = EXPERIMENTALLY OBSERVED NA2S04 DEPOSTION RATE (MG/HR) 

XJ(I) = MOLE FRACTION OF SPECIES I IN FLAME 

XI = TURINXRE/100,F(TURB) INTERMEDIATE 

XWCI) = MOLE FRACTION OF SPECIES I AT WALL 

X ( 4 ) = MOLE FRACTION OF DOMINANT SPECIES IN MIXTURE 

YNAJ = TOTAL SODIUM MOLE FRACTION IN FLAME 

YNAW = TOTAL SODIUM MOLE FRACTION AT STATION W 

YSJ = TOTAL SULFUR MOLE FRACTION IN FLAME 

YSW = TOTAL SULFUR MOLE FRACTION AT STATION W 


CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

c 

c 
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*0016300 

0016400 

O016500 

0016600 

0016700 

0016800 

0016900 

0017000 

0017100 

0017200 

0017300 

0017400 

0017500 

0017600 

0017700 

0017800 

0017900 

0018000 

0018100 

0013200 

0013300 

0018400 

0018500 

0018600 

0013700 

0018800 

0018900 

0019000 

0019100 

0019200 

0019300 

0019400 

0019500 

0019600 

0019700 

0019300 

0019900 

0020000 

0020100 

0020200 

0020300 

0020400 

0020500 

0020600 

0020700 

0020800 

0020900 

0021000 

0021100 

0021200 

0021300 

0021400 

0021500 

0021600 


C SODIUM VERSION OF CHEMICALLY FROZEN BOUNDARY LAYER THEORY 
C USED FOR SODIUM SULFATE DEPOSITION RATE CALCULATIONS 
C 

C MAIN PROGRAM 
C 

REAL LAMIX,LAMDA,LEW,LEJ,NU,LW 
REAL M(7) 

DIMENSION V(7),XWC7),XJ(7),SM(7),ALPHA(7),CMHC7) 

CONMON/VA/DATOS ( 3,5) 

LOGICAL FSORET 
COMMON FSORET 
C0NM0N/GFR/X(4) 

DATA V/ 1.0 , 1.0 , 2.0 , 0.0 , 0.0 , 1.0 , 0.0 / 

C 

WRITEC6 ,490 ) 

490 FORMAT ( 10X, f CFBL THEORY FOR SODIUM SULFATE DEPOSITION RATES/) 
WRITEC6, 900 ) 

900 FORMAT ( 1 OX, f ALL GAS PROPERTIES PERTAIN TO TJ,PJ f ,/) 

C 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
C 

C PAY SPECIAL ATTENTION TO THE FORMAT OF THE INPUT PARAMETERS 

C SEE 'VARIABLE NAME LIST* ABOVE FOR DESCRIPTION OF VARIABLES 

C 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

c 

C DOMINANT SPECIES IN THE MIXTURE 

C l=NITROGEN , 2=0XYGEN , 3=WATER , 4=CARB0N DIOXIDE 
C 

C TRACE SPECIES IN THE MIXTURE 

C l=NAOH , 2 = NA , 3=NA2S04 , 4=SQ2 , 5=S03 , 6=NACL , 7=H2S 

C 

C 

C INPUT : 1 
C 

READ(5,1)IRUN, ITYPE, F,WA,TW,TO,PO, PJ, DC, WOBS, FSORET 

1 FORMAT (2I2,8F9.0,3X,L1) 

C 

C INPUT : 2 
C 

READC 5 , 998 ) DI AW, LW,DJ,TURB, TURIN, TURL, SHAPE, DIV 
998 FORMAT (8F1 0 . 0 ) 

C 

C INPUT : 3 
C 

DO 502 1=1,7 

READC 5 , 2 ) XW ( I ) , X J ( I ) , CMH Cl) 

2 FORMAT (2E15.4,F10.0) 

502 CONTINUE 

C 

C OUTPUT OF INPUT PARAMETERS 
C 

WRITEC 6 , 978 ) IRUN , ITYPE, F, TW, TO , PO , PJ , DC, WOBS , FSORET 
978 FORMAT (10X, 1 IRUN = f , 12 ,/ , 10X, 'ITYPE = ' , 12,/, 10X, ' F =', 
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0021700 
0021800 
0021900 
0022000 
0022100 
0022200 
0022300 
0022400 
0022500 
0022600 
0022700 C 
0022800 C 
0022900 C 
0023000 
0023100 
0023200 
0023300 
0023400 
0023500 
0023600 
0023700 
0023800 
0023900 
0024000 
0024100 
0024200 C 
0024300 C 
0024400 C 
0024500 C 
0024600 
0024700 
0024800 
0024900 
0025000 
0025100 
0025200 
0025300 C 
0025400 C 
0025500 C 
0025600 C 
0025700 C 
0025800 C 
0025900 
0026000 
0026100 
0026200 
0026300 
0026400 
0026500 C 
0026600 C 
0026700 C 
0026800 
0026900 
0027000 


*F10.6,/,10X, T TW (K) = f ,F10.3,/,10X, 

#'TO (K) =',F10.3,/,10X,'PO (ATM) = ', F10 . 6 ,/, 10X, * PJ (ATM) =»,F10.6 
*,/,10X, f DC = ' ,F9.6/10X, 'WOBS (MG/HR) = ' , F10 . 4 ,/, 10X, f FSORET =',L1) 
UIR I T E ( 6 , 988 ) WA , DI AW, LW, DJ , TURIN , TURL , SHAPE, DIV 
988 FORMAT ( 10X, f WA (G/SEC) = * , F10 . 5 , / , 1 OX, 1 DI AW (CM) = ? ,F10.5, 
*/,10X,’LW (CM) = f ,F10.5,/,10X, T DJ (CM) = f , F10 . 5 ,/, 10X, 

* f TURIN =’,F8.4, f PERCENT', /, 10X, 'TURL (CM) = ' , F10 .5,/, 10X, 

* 1 SHAPE = ' ,F7.4,/,10X, f DIV =',F7.4) 

WRITE(6,939) 

939 FORMAT(/,10X, *I=1=NA0H , I=2=NA , I=3=NA2S04 , I=6=NACL',/) 

ASSIGN E/K, SIGMA AND MOLECULAR WEIGHT OF DOMINANT SPECIES 

DATOS ( 1 , 1 ) =71 . 4 
DAT0S(2,1)=3.798 
DATOS(3,1)=28.0134 
DATOS(1,2)=106.7 
DAT0S(2,2)=3.467 
DAT0S(3,2)=31. 9988 
DATOS ( 1,3) =356 . 0 
DAT0S(2,3)=2.649 
DAT0S(3,3)=18. 0152 
DAT0S(1,4)=195.2 
DAT0SC2,4)=3. 941 
DATOS(3,4)=44.0098 

THERMAL DIFFUSION FACTOR OF TRACE SPECIES 

I==NAOH, I=2=NA, I=3=NA2S04,I=4=S02,I=5=S03,I=6=NACL,I=7=H2S 


ALPHA C 1 ) =0 . 1592*(l.-429.5/TW) 

ALPHA(2)=-0. 0607* (1.-551. 2/TW) 

ALPHA ( 3 )=0 .7828* (1.-476 .2/TW) 

ALPHA (4) =0 .3131*(1 .-198.2/TW) 

ALPHA ( 5) =0. 3877* ( 1 . -184 . 3/TW) 

ALPHA(6)=0.3462*(1 .-420.0/TW) 

ALPHA(7)=0 .0592*(l.-79.38/TW) 

COMPUTE THE AVERAGE MOLECULAR WEIGHT OF MIXTURE 
(NEGLECT THE TRACE SPECIES) 

1= NITROGEN* 2 = OXYGEN, 3=WATER, 4=CARB0N DIOXIDE 
MOLE FRACTIONS OF DOMINANT SPECIES AFTER COMBUSTION 

DEN=1.+1.0331*F 
X(1)=0. 7905/DEN 
X(2)=(0. 2095-3. 0993*F)/DEN 
X(3)=2. 0662* F/ DEN 
X(4)=X(3) 

AM=DAT0S(3,1)*X(1)+DAT0S(3,2)*X(2)+DAT0S(3,3)*X(3)+DAT0S(3,4)*X(4) 

CALCULATE JET EXIT TEMPERATURE, (K) , AND VELOCITY, (CM/SEC) 

CALL CALMIX(TO,AM,CMIX) 

GM=CMIX/ ( CMIX-1 . 9872/AM) 

GMR= ( GM-1 . )/GM 
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P 027 1 0 0 
0027200 
0027220 
0027300 
0027400 C 
0027500 C 
0027600 C 
0027700 
0027800 
0027900 
0028000 
0028100 
0028200 
0028300 
0028400 
0028500 
0028600 C 
0028700 C 
0028800 C 
0028900 
0029000 
0029100 
0029200 
0029300 C 
0029400 C 
0029500 C 
0029600 
0029700 
0029800 
0029900 
0030000 
0030100 
0030200 
0030300 
0030400 
0030500 
0030600 
0030700 
0030800 
0030900 
0031000 
0031100 
0031200 
0031300 
0031400 
0031500 
0031600 
0031700 
0031800 
0031900 
0032000 
0032100 
0032200 
0032300 


AJ=3. 14159* DJ*DJ/4. 

TJ=TO*(PJ/PO)**GMR 
IFCDC.EQ.0.00) DC = 1 .00 

UJ=( (82 . 057*TJ )/( PJ*AM) )*(WA*(1*+F) )/( AJ*DC) 

CALCULATE MIXTURE PROPERTIES AT STATIONS J AND W 

RHO J=PJ*AM/T J/82 . 057 
RHOW=PJ*AM/TU/82.057 
CALL CALMIX(TJ,AM,CMIX) 

CALL CALMIX(TW,AM,CP) 

CALL MIXPRO(TJ,ETAMIX,LAMIX) 

CALL MIXPROCTW,ETA,LAMDA) 

APHMIX=L AMIX/CMIX/RHQJ 

APHW=LAMDA/CP/RHOW 

PR=CNIX*ETAMIX/LAMIX 

CALCULATE REYNOLDS NUMBER 

IFC SHAPE . EQ . 0 . 0 ) SHAPE=1.0 
IFCDI V . EQ . 0 . 0 ) DIV=1.0 
REJ=RHO J*UJ*DI AW/ETAMIX 
RE=REJ*SHAPE*DIV 

GIVE VALUES OF E/K, SIGMA AND MOLEC. WEIGHT FOR TRACE SPECIES I 

DO 10 1=1,7 
25 I F( I . GT . 1 ) GO TO 5 
DATOS(1,5)=1962.0 
DATOS(2,5)=3.804 
DAT0S(3,5)=39.9971 
GO TO 11 

5 IFCI.GT.2) GO TO 6 
DATOSC 1,5) =1375 . 

DAT0S(2,5)=3.567 
DAT0SC3,5)=22. 98977 
GO TO 11 

6 IFCI.GT.3) GO TO 7 
DATOSC 1,5) =2221 . 

DATOSC 2 , 5 ) =5 . 0 0 
DATOSC3,5)=142.0371 
GO TO 11 

7 IFCI.GT.4) GO TO 8 
DAT0SC1,5)=335.4 
DATOSC 2 , 5 ) =4 . 1 12 
DATOS C 3 , 5 ) =64 . 0588 
GO TO 11 

8 IFCI.GT.5) GO TO 9 
DATOSC 1,5) =431 . 4 
DATOS(2,5)=4.207 
DA T0SC3, 5) =80.0582 
GO TO 11 

9 IFCI.GT.6) GO TO 12 
DAT0S(1,5)=1989. 
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0032400 DAT0S(2,5)=4.186 

0032500 DAT0S(3,5)=58.4428 

0032600 GO TO 11 

0032700 12 CONTINUE 

0032800 DAT0S(1,5)=301.1 

0032900 DAT0SC2,5)=3.623 

0033000 DAT0SC3, 55=34.0758 

0033100 11 CONTINUE 

0033200 C 

0033300 C DIFFUSION COEFFICIENT OF TRACE SPECIES AT STATIONS J AND M 
0033400 C 

0033500 CALL DIFC T J , P J , DA ) 

0033600 CALL DI FC TW, P J , DW> 

0033700 C 

0033S00 C LEWIS NUMBER OF TRACE SPECIES AT STATIONS J AND U 
0033900 C 

0034000 LEJ=DA/APHMIX 

0034100 LEW=DW/APHW 

0034200 C 

0034300 C SCHMIDT NUMBER OF TRACE SPECIES AT STATION J 
0034400 C 

0034500 SC=ETAMIX/RHO J/DA 

0034600 C 

0034700 C MASS TRANS. NUSSELT NO . (PERIMETER-AVERAGED, CYLINDERICAL COLLECTOR) 
0034800 C 

0034900 CALL NUM( TW, TO ,’RE, SC, NU) 

0035000 C 

0035100 C COMPUTE F(SORET) 

0035200 C 

0035300 CALL TMDFC I , TW, TO , LEW, ALPHA, TAU , THERM) 

C035400 C 

0035500 C COMPUTE F(TURB) 

0035600 C 

0035700 CALL TURBL ( RE, TURB , TURIN, TURL , DIAW) 

0035800 C 

0035900 C CALCULATE NA2S04 DEPOSITION RATE CONTRIBUTED FROM SPECIES I 
0036000 C BY CFBL FORMULATION. 

0036100 C 

0036200 IF(CMHCI).EQ.O.O) CMH(I)=1.0 

0036300 SM(I)=TURBXDA*RHOJXTHERMX((XJ(I)-XW(I))+XW(I)XTAU/THERM# 

0036400 x(LEW/LEJ)*xo • 6*CMIX/CP/CMH(I) )*NU*DATOS(3,5)/AM/DIAW 

0036500 IFCV(I).EQ.O.O) GO TO 10 

0036600 AW=3. 14159*DIAW*LW 

0036700 M(I)= A W*3. 600 E6 *0.5*142. 037 1/DAT0S(3,5)*SM(I)XV(I) 

0036800 C 

0036900 C PRINT INTERIM OUTPUT 
0037000 C 

0037100 IF(I.GT.l) GO TO 546 

0037200 IF(FSORET) GO TO 100 

0037300 WRITE(6,544) 

0037400 544 FORMAT(/, 10X, 'CALCULATIONS DO NOT INCLUDE THERMAL DIFFUSION',/) 

0037500 100 WRITE(6,545) 

0037600 545 FORMAT ( T2 , • I ' , Til , ' D( I ) ' ,T25, 'NUC I ) ' , T40 , *SC( I ) ' , T55, 'X(I) ,- 

0037700 CJ’ ,T70, » X( I ) , W f , T85 , ' TAUC I ) ' , T100 , ' F( SORET ) , I ' , T116 , 'M( I ) ’ ) 
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0037800 

p 0 37 90 0 

0038000 

0038100 C 

0038200 C 

0038300 C 

0 0 3S40 0 C 

0038500 

0038600 

0033700 

0038800 

0038900 C 

0039000 C 

0039100 C 

0039200 

0039300 C 

0039400 C 

0039500 C 

0039600 

0039700 

C 0 3980 0 

0039900 

0040000 

0040100 

0040200 

0040300 

0040400 C 

0040500 C 

0040600 C 

0040700 

0040300 C 

0040900 C 

00^1000 C 

0041100 

0041200 

0041300 C 

0041400 C 

0041500 C 

0041600 

0041700 

0041800 

0041900 

0042000 

0042100 

0042200 

0042300 

0042400 

0042500 

0042600 

0042700 

0042800 

0042900 

0043000 

0043100 


546 WRITEC6 > 550 ) I,DA,NU,SC,XJ(I),XWCI),TAU,THERM,MCI) 

55 0 FORMAT (IX, II, 6 El 5.5* FI 5. 5, 2 El 5. 5*/) 

10 CONTINUE 

COMPUTE TOTAL SODIUMCNA) AND SULFUR(S) MASS FLUX 
SUM=J(NA), SMS=J(S) 

SUM=SMC1)*22. 98977/39. 9971+SMC2)*! . 0+SMC3)*2.*22. 98977 
// 142. 0371+ SM (6)^22.98977/58. 44 28 
SMS=SM( 3 )*32. 06/142. 0371+SMC4 )*32. 06/64 .0588+SM(5)K 
*32. 06/80. 0582+SMC 7) *32. 06/34. 0758 

COMPUTE SODIUM TO SULFUR MOLAR FLUX RATIO 

RATI0=SUM/SMS*32 .06/22.98977 

COMPUTE TOTAL NA AND S MOLE FRACTIONS AT STATIONS J AND U 

YNAJ=XJ(1)+XJ(2)+(XJ(3)*2)+XJ(6) 

YNAU=XW(1)+XW(2)+(XW(3)*2)+XW(6) 

YS J =XJ ( 3 ) +XJ ( 4 )+XJ ( 5 )+XJ (7 ) 

YSW=XW(3)+XW(4)+XW(5)+XW(7) 

WRITE (6,211) 

211 FORMAT (//,T8, 'YNAJ * ,T24, f YNAW f ,T40, f YSJ f ,T57, f YSW f ) 

WRITE(6, 212) YNAJ, YNAW,YSJ,YSW 

212 F0RMATC4E16.6,/) 

CALCULATE NA2S04 DEPOSITION RATE (MG/HR) 

WPRED=M(1)+M(2)+M(3)+M(6) 

COMPUTE PERCENT ERROR BETWEEN M(PRED) AND M(OBS) 

ERR=WPRED-WOBS 
ERROR = ERR/ WOBS*I 00 . 

PRINT RESULTS 

WRITE(6,4)AM,TURB,RE,SUM,SMS,RATI0,WPRED,ERR0R,GM,PR,X(1),XC2),X(3), 

*X(4),TJ 

4 FORMAT (/ , 1 OX, 1 AM = f , FI 0 . 5,/ , 1 OX, f TURB =* , F10 . 5,/, 10X, *RE =*,F10.3,/, 
*1 OX, 1 SUM (G/CMXX2/SEC) = f , EI3 . 5 ,/, 10X, 1 SMS (G/CM**2/SEC) = f ,E13.5,/, 
*1 OX , f RATIO =' ,F10.4,/,10X, 

* 1 PREDICTED NA2S04 DEPOSITION RATE (MG/HR) = 1 , F10 . 5,/, 10X, 

* f ERROR (X) =' ,F10.4,/,10X, f GM = ' , F10 .5,/, 10X, 1 PR = 1 , F10 . 6 ,/, 10X, 
*'X(N2) = f ,F10.6,/,10X, *X( 02) = f , F10 . 6 ,/, 1QX, 1 XCH20 ) = ' , F10 .6 ,/, 10X, 
#'X(C02) = 1 ,F10.6,/,10X, f TJ (K) = f ,F10.4,/) 

WRITEC6 , 950 ) 

950 FORMAT (/,T5, ' RHO J ( G/CM**3 ) 1 >T21, ' UJ ( CM/S) ' , T37 , f ETAMIX(POISE) 

C* ,T52, r LAMIX(CAL/CM/K/S) f ,T73, T CMIXCCAL/G/S) 1 ) 

WRITE(6,951)RH0J,UJ,ETAMIX,LAMIX,CMIX 

951 FORMAT (4E16.6,4X,E16.6,//) 

STOP 

END 
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G 04320 0 
0043300 
0043400 
0043500 
0043600 
0043700 
0043300 
0043900 
0044000 
0044100 
0044200 
0044300 
0044400 
0044500 
0044600 
0044700 
0044800 
0044900 
0045000 
0045100 
0045200 
0045300 
0045400 
0045500 
0045600 
0045700 
0045800 
0045900 
0046000 
0046100 
0046200 
0046300 
0046400 
0046500 
0046600 
0046700 
0046800 
0046900 
0047000 
0047100 
0047200 
0047300 
0047400 
0047500 
0047600 
0047700 
0047S00 
0047900 
0048000 
0048100 
0048200 
0048300 
0048400 
0048500 


C 

cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx* 

c 

SUBROUTINE PARAM ( J , EPSOVK, SIG2, M) 

c 

C SUBROUTINE PARAM GIVES THE VALUES OF E/K, SIGMA AND MOLECULAR WEIGHT 
C 

REAL M 

C0MM0N/VA/DAT0S(3,5) 

EPSOVK=DATOS(l,J) 

SIG2=DAT0S(2, J) 

M=DATOS ( 3 , J ) 

RETURN 

END 

C 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

c 


c 

c 

c 

c 


0.16022128E-02 , -0 . 62936893E-06 
0 . 11441022E-Q9 ,-0 . 78057465E-14 
0 . 36 122139E+01 
0.33749008E-10 
0 . 26340654E+01 , 0 . 31121899E-02 , -0 . 90278449E-06 
0. 12673054 E- 09 ,-0 . 69164732E-14 , 

Q.3Q981719E-02 


0 . 74853166 E-03 , -0 . 19820647E-06 
,-0. 23907374 E- 14 


0 . 44608041E+01 


SUBROUTINE CAL ESL C J , T , CJ ) 

THIS SUBROUTINE CALCULATES DIMENSIONLESS HEAT CAPACITY FROM A 
CURVE FIT VALID FOR 300K < T < 1000K AND 1000K < T < 5000K 

DIMENSION C0EF1(5,4),C0EF2C5,4) 

DATA COEF1/ 0 . 28532899E+01 
x 
x 

X 
X 
X 
X 

* 0.22741325E-09 , -0 . 15525954E-13 

DATA C0EF2/ 0 . 37 044177 E+0 1 , -0 . 14218753 E-02 

* -0.12028885E-08 

* 0 . 37837135E+01 

* -0 . 98189101E-08 

* 0 . 41675564E+01 

* -0 .48670871E-08 

* 0 .24007797E+01 

x 0.20021861E-08 

IFCT.GT.1000. ) GO TO 10 
CJ=C0EF2(5, J)XT**4+C0EF2(4, J)XT**3+C0EF2(3, J)*T*X2+C0EF2(2, J)*T 
*+C0EF2(l,J) 

GO TO 20 

10 CJ=C0EF1C5, J)XTX*4+C0EF1(4, J)*T*X3+C0EF1(3, J)XT**2+C0EF1(2, J)*T 
x + COEFHl, J) 

20 RETURN 
END 


, “0 . 13954677 E- 13 
/ “0, 30233634 E“ 02 
, 0 . 330 3 1825E-1 1 
,-0. 18106868 E- 02 
, 0.15284144E-11 
, 0.87350 957E-02 
, 0.63274039E-15 


,-0.1239257 IE- 05 
/ 

, 0.28670392E-05 
l 0 . 99492751 E-05 
l 0 . 59450878E-05 
^-0. 66070878 E-05 


CXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

c 


SUBROUTINE CALMIXC T , AM, CMIX) 

C 

C SUBROUTINE CALMIX PROVIDES MIXTURE HEAT CAPACITY IN CAL/(G*K) 
C 
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0 0436 0 0 
0048700 
0043800 
0 04S 90 0 
0049000 
0049100 
0049200 
0049300 
0049400 
0049500 
0049600 
0049700 
0049800 
0049900 
0050000 
0050100 
0050200 
0050300 
0050400 
0050500 
0050600 
0050700 
0050800 
0050900 
0051000 
0051100 
0051200 
0051300 
0051400 
0051500 
0051600 
0051700 
0051800 
0051900 
0052000 
0052100 
0052200 
0052300 
0052400 
0052500 
0052600 
0052700 
0052800 
0052900 
0053000 
0053100 
0053200 
0053300 
0053400 
0053500 
0053600 
0053700 
0053800 
0053900 


COMMON/GFR/XC 4 ) 

C = 0. 

DO 8 J=l,4 

CALL CALESL(J,T,CJ) 

C=C+CJ*X(J)*1.9872 
8 CONTINUE 
CMIX=C/AM 
RETURN 
END 
C 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

C 

SUBROUTINE COLINTCII, JJ,T,OVDT,ETA,CONDT,OMEDT) 

C 

C THIS SUBROUTINE GIVES COLLISION INTEGRALS , VISCOSITY AND THERMAL 
C CONDUCTIVITY OF SPECIES USING CHAPMAN-ENSKOG THEORY 
C 

REAL M 

DIMENSION 0ME(79),ALPA(79),0MED(79) 


DATA 

ALPA/ 0.300 , 0.350 , 

0.400 , 

0.450 , 0.500 , 0. 

550 , 0.600 

9 - 

X 

0.650 , 0.700 , 

0.750 , 

0.800 , 0.850 > 0. 

900 , 0.950 


X 

1.000 , 1.050 , 

1.100 , 

1.150 , 1.200 , 1. 

250 , 1.300 

f - 

X 

1.350 , 1.400 , 

1.450 , 

1.500 , 1.550 , 1. 

600 , 1.650 

, - 

X 

1.700 , 1.750 , 

1.800 , 

1.850 , 1.900 , 1. 

950 , 2.000 

f - 

X 

2.100 , 2.200 , 

2.300 , 

2.400 , 2.500 , 2. 

600 , 2.700 

f - 

X 

2.800 , 2.900 , 

3.000 , 

3.100 , 3.200 , 3. 

300 , 3.400 

, - 

X 

3.500 , 3.600 , 

3.700 , 

3.800 , 3.900 , 4. 

000 , 4.100 

> - 

X 

4.200 , 4.300 , 

4.400 , 

4.500 , 4.600 , 4. 

700 , 4.800 

9 - 

X 

4.900 , 5.000 , 

6.000 , 

7.000 , 8.000 , 9. 

000 , 10.00 

9 - 

X 

20.00 , 30.00 , 

40.00 , 

50.00 , 60.00 , 70 

.00 , 80.00 

9 - 

X 

90.00 , 100.0 / 





DATA 

OME/ 2.7850 , 2.6280 

, 2.4920 

, 2.3680 , 2.2570 

, 2.1560 , 

- 

X 

2.0650 , 1.9820 

, 1.9080 

, 1.8410 , 1.7800 

, 1.7250 , 

- 

X 

1.6750 , 1.6290 

, 1.5870 

, 1.5490 , 1.5140 

, 1.4820 , 

- 

X 

1.4520 , 1.4240 

, 1.3990 

, 1.3750 , 1.3530 

, 1.3330 , 

- 

X 

1.3140 , 1.2960 

, 1.2790 

, 1.2640 , 1.2480 

, 1.2340 , 

- 

X 

1.2210 , 1.2090 

, 1.1970 

, 1.1860 , 1.1750 

, 1.1560 , 

- 

X 

1.1380 , 1.1220 

, 1.1070 

, 1.0930 , 1.0810 

, 1.0690 , 

- 

X 

1.0580 , 1.0480 

, 1.0390 

, 1.0300 , 1.0220 

, 1.0140 , 

- 

X 

1.0070 , 0.9999 

, 0.9932 

, 0.9870 , 0.9811 

, 0.9755 , 

- 

X 

0.9700 , 0.9649 

, 0.9600 

, 0.9553 , 0.9507 

, 0.9464 , 

- 

X 

0.9422 , 0.9382 

, 0.9343 

, 0.9305 , 0.9269 

, 0.8963 , 

- 

X 

0.8727 , 0.8538 

, 0.8379 

, 0.8242 , 0.7432 

, 0.7005 , 

- 

X 

0.6718 , 0.6504 

, 0.6335 

, 0.6194 , 0.6076 

, 0.5973 , 

- 

X 

0.5882 / 





DATA 

OMED/ 2.6620 , 2.4760 

, 2.3180 

, 2.1840 , 2.0660 

, 1.9660 , 

- 

X 

1.8770 , 1.7980 

, 1.7290 

, 1.6670 , 1.6120 

, 1.5620 , 

- 

X 

1.5170 , 1.4760 

, 1.4390 

, 1.4060 , 1.3750 

, 1.3460 , 

- 

X 

1.3200 , 1.2960 

, 1.2730 

, 1.2530 , 1.2330 

, 1.2150 , 

- 

X 

1.1980 , 1.1820 

, 1.1670 

, 1.1530 , 1.1400 

, 1.1280 , 

- 

X 

1.1160 , 1.1050 

, 1.0940 

, 1.0840 , 1.0750 

, 1.0570 , 

- 

X 

1.0410 , 1.0260 

, 1.0120 

, 0.9996 , 0.9878 

, 0.9770 , 

- 

X 

0.9672 , 0.9576 

, 0.9490 

, 0.9406 , 0.9328 

, 0.9256 , 

- 

X 

0.9186 , 0.9120 

, 0.9058 

, 0.8998 , 0.8942 

t 0 . 8888 p 
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0054000 
0054100 
0054200 
0054300 
0054400 
0054500 
0054600 
0054700 
0054800 
0054900 
0055000 
0055100 
0055200 
0055300 
0055400 
0055500 
0055600 
0055700 
0055800 
0055900 
0056000 
0056100 
0056200 
0056300 
0056400 
0056500 
0056600 
0056700 
0056800 
0056900 
0057000 
0057100 
0057200 
0057300 
0057400 
0057500 
0057600 
0057700 
0057800 
0057900 
0058000 
0058100 
0058200 
0058300 
0058400 
0058500 
0058600 
0058700 
0058800 
0058900 
0 0 5 9 0 0 0 
0059100 
0059200 
0059300 


0.8836 , 

0.8788 , 

0.8740 , 

0.8694 , 

0.8652 , 

0.8610 , - 


0.8568 , 

0.8530 , 

0.8492 , 

0.8456 , 

0.8422 , 

0.8124 , - 

* 

0.7896 , 

0.7712 , 

0.7556 , 

0.7424 , 

0.6640 , 

0.6232 , - 


0.5960 , 

0.5756 , 

0.5596 , 

0.5464 , 

0.5352 , 

0.5256 , - 



* 0.5170 / 

I F( J J . EQ . 0 ) GO TO 67 

CALL PARAM(JJ,EOKJ,SIGL,AJJ) 

CALL PARAM(II,EOKI,SIGK,AII) 

M=(AII+AJJ)/2. 

SIG2=(SIGL+SIGK)/2. 

EOVKIJ=SQRT(EOKJXEOKI) 

ADT=T/EOVKIJ 
GO TO 68 

67 CALL PARAM(II,EPS0VK,SIG2,M) 

ADT=T/EPSOVK 

68 J - 3 

69 IF(ADT-ALPA(J))112,111,110 

110 J = J + 1 

IF( J . LE.79)GO TO 69 
GO TO 204 

111 OVDT-OME(J) 

OMEDT = OMED( J ) 

GO TO 200 

112 B=(ADT-ALPA(J-1))/(ALPA(J)-ALPA(J-1)) 

OVDT=B*(B-l. )*OME( J-2)/2.+Cl.-BX*2)X0MECJ-l)+BX(B+l. )*OMECJ)/2. 
OMEDT =BX(B-1 . )*OMED( J-2)/2.+(l .-BX*2)X0MED( J-1)+B*(B+1. )*OMED( J)/2- 
x . 

GO TO 200 

204 OVDT=0.5882X((ADT/100.)XX(-0.17)) 

OMEDT =0.5170X( ( ADT/1 0 0 . )*X(-0 . 15) ) 

200 ETA=26 .693*SQRT(MXT)/( (SIG2)**2X0VDT)X1. E-06 
CALL CALESL(II,T,CR) 

CONDT-Cl. 9872/M) *C 15. /4.+1. 32* ( CR-5. /2. ) )*ETA 
RETURN 
END 
C 

C 

SUBROUTINE MIXPRO ( T , ETAMIX, L AMIX) 

C 

C THIS SUBROUTINE PROVIDES MIXTURE VISCOSITY AND THERMAL CONDUCTIVITY 
C 

REAL LAMP , L AMIX 
DIMENSION VISCC 4 ) > LAMP (4) 

COMMON/VA/DATOS ( 3,5) 

C0MM0N/GFR/XC4) 

ETAMIX=0. 

L AMIX-0 . 

11 = 0 

DO 40 J = l,4 
N = J 

CALL COLINTCN,Il,T,OVDT,ETAMJ,CQNDT,OMEDT) 

VISCC J)=ETAMJ 
LAMPC J)=CONDT 
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0059400 

0059500 

0059600 

0059700 

0059800 

0059900 

0060000 

0060100 

0060200 

0060300 

0060400 

0060500 

0060600 

0060700 

0060800 

0060900 

0061000 

0061100 

0061200 

0061300 

0061400 

0061500 

0061600 

0061700 

0061800 

0061900 

0062000 

0062100 

0062200 

0062300 

0062400 

0062500 

0062600 

0062700 

0062800 

0062900 

0063000 

0063100 

0063200 

0063300 

0063400 

0063500 

0063600 

0063700 

0063800 

0063900 

0064000 

0064100 

0064200 

0064300 

0064400 

0064500 

0064600 

0064700 


40 CONTINUE 
DO 18 1=1,4 
PROV=0 . 

DO 28 J=1 , 4 

G=1./(2.*SQRT(2.))X((1.+DATOS(3,I)/DATOSC3,J))XX(-0.5))XC1. 
x +SQRT(VISC(I)/VISCCJ) )XSQRT(SQRT(DAT0S(3, J)/DAT0S(3,I))))XK2 
PROV=PROV+GXX( J) 

28 CONTINUE 

ETAMIX=ETAMIX+XCI)XVISC(I)/PRQV 
L AMIX=L AMIX+XC I ) XL AMP ( I )/PROV 
18 CONTINUE 
RETURN 
END 
C 

CXHXKXKKK*KKXXXHHXXXXKKKKKK#XXXXKKXXKXXKKHKXKXXXHXXXKXXXXKXXXKXXXXXKXXXXXXKKXXKX*HX 

C 


SUBROUTINE COEDIFCT, P , I , J , D) 

C 

C THIS SUBROUTINE CALCULATES BINARY DIFFUSION COEFFICIENT 
C 

CALL COLINTCI, J , T, OVDT, ETA , CONDT , OMEDT) 

CALL PARAMCI,EPSOVK,SIGI,WI) 

CALL PARAMC J , EPSOVK, SIGJ, WJ) 

SIG=1 ./2.x(SIGI+SIGJ) 

D-0 . 0018583XSQRT C (TX*3)X(1 ,/WI+l ,/WJ ) )/CPX0MEDTX(SIGXX2) ) 

RET URN 
END 
C 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

SUBROUTINE DIF(T,P,DA) 

c 

C IN B GAS TIN X DIF GIVES DIFFUSI0N COEFFICIENT OF TRACE SPECIES K 

C 

C0NM0N/GFR/XC4) 

DEN=0 . 0 
DO 101 J=l,4 
K = 5 

CALL COEDIFCT, P,K,J,DIFF) 

DEN=DEN+X(J)/DIFF 
101 CONTINUE 
DA=1./DEN 
RETURN 
END 
C 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

Vf 


c 

c 

c 

c 


SUBROUTINE TMDFC I , TW, TO , LEW, ALPHA , TAU, THERM) 


SUBROUTINE TMDF CALCULATES THERMAL 
LABELLED AS 'THERM'. IF FSORET IS 


DIFFUSION EFFECT, 
FALSE, THERM=1 . 0 


F(SORET), 


LOGICAL FSORET 


29 



0064800 COMMON FSORET 

0064900 REAL LEW 

0065000 DIMENSION ALPHAC7) 

0065100 IF(FSORET) GO TO 10 

0065200 THERM=1 . 0 

0065300 T AU = 0 , 0 

0065400 GO TO 20 

0065500 10 TAU=ALPHA(I)KLEW**0.4H(TO-TW)/TW 

0065600 ET=EXP(-TAU) 

0065700 THERM=TAU/( 1 . -ET) 

0065800 20 RETURN 

0065900 END 

0066000 C 

0 066 ICO CX*XKKXX*KKKXXXXMXXXXX*XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXX 
0066200 C 

0066300 SUBROUTINE NUM(TW,TO,RE,SC,NU) 

0066400 C 

0066500 C SUBROUTINE NUM CALCULATES PERIMETER AVERAGED NUSSELT NUMBER 
0066600 C INCLUDING EFFECTS OF MACH NUMBER AND VARIABLE PROPERTIES 
0066700 C 

0066800 REAL NU 

0066900 NU=(0 . 4 0 *SQRT ( RE) +0 . 06*RE*X(2 . /3 . ) )XCSCXX0 .4)* 

0067000 1(TO/TW)**(-0.04) 

0067100 RETURN 

0067200 END 

0067300 C 

0067400 c xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
0067500 C 

0067600 SUBROUTINE TURBL ( RE, TURB , TURIN , TURL , DI AW) 

0067700 C 

0067800 C SUBROUTINE TURBL CALCULATES EFFECT OF MAINSTREAM TURBULENCE 

0067900 C 

0068000 REAL LOD 

0068100 I FC TURIN . NE . 0 . 0 ) GO TO 5 

0068200 IFCTURB.EQ.O.O) GO TO 10 

006S300 GO TO 30 

0068400 5 LOD'TURL/DI AW 

0068500 IFCLOD.GT.2.0) PSI =0 . 124* C LOD-1 1 . 0 ) XX2 . +2 . 

0068600 IFCLOD.LE.2.0) PSI =***4 •0*(LOD-1.75)**2.+12.25 

0068700 PSI=PSI/1000. 

0068800 XI Z TURINXRE/100. 

0068900 IFCXI-100. )1,2,2 

0069000 1 PHI=12.375x(l.-(l.-XI/100. )XX1.5) 

0069100 GO TO 20 

0069200 2 PH 1=9. 0+0.03375^X1 

0069300 20 TURB=1.+PSI*PHI 

0069400 GO TO 30 

0069500 10 TURB=1.0 

0069600 30 RETURN 

0069700 END 

0069800 C 


0069900 CXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX END OF PROGRAM XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
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